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Abstract 

In this paper we obtain an analytical solution of the relativistic Boltzmann equation under 
the relaxation time approximation that describes the out-of-equilibrium dynamics of a radially 
expanding massless gas. This solution is found by mapping this expanding system in flat spacetime 
to a static flow in the curved spacetime AdS 2 <8> S 2 • We further derive explicit analytic expressions 
for the momentum dependence of the single particle distribution function as well as for the spatial 
dependence of its moments. We find that this dissipative system has the ability to flow as a perfect 
fluid even though its entropy density does not match the equilibrium form. The non-equilibrium 
contribution to the entropy density is shown to be due to higher order scalar moments (which 
possess no hydrodynamical interpretation) of the Boltzmann equation that can remain out of 
equilibrium but do not couple to the energy-momentum tensor of the system. Thus, in this system 
the slowly moving hydrodynamic degrees of freedom can exhibit true perfect fluidity while being 
totally decoupled from the fast moving, non-hydrodynamical microscopic degrees of freedom that 
lead to entropy production. 
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I. INTRODUCTION 


The Boltzmann equation describes the underlying microscopic dynamics of dilute classical 
gases [l|. It is widely employed to model a variety of nonequilibrium phenomena in several 


areas of physics such as the dynamics o 
of ultrarelativistic heavy ion collisions 


■ the hot hadronic matter produced in the late stages 
, y] , some aspects of the expansion of our universe 
in cosmology applications jd], the description of micro and nano-flows j^, among others. 

In addition to these app lications, exact solutions of the relativistic generalization of the 


Bolt...... equation flQ i.. tl.e aeWtio.. ti..e app.oxi..atio.. gff 


have been recently 


employed to improve our understanding of the domain of applicabilit y o f relativistic dissi¬ 
pative fluid dynamics in the context of relativistic heavy ion collisions jl0l-ll4l|. Even though 
less complete, the Anderson-Witting-Boltzmann (AWB) equation and its solutions can be 
used to understand certain properties of solutions of the Boltzmann equation itself, as well 
as its hydrodynamic limit. 

Analytic solutions of the relativistic Boltzmann equation are extremely rare (see 151 
for the first analytical solution in an expanding background). The same can be said even 
for simplihed versions of the relativistic Boltzmann__equation, such as the AWB equation. 

was derived for a conformal system 


, (for an extension 


Recently, an exact solution of the AWB equation 
undergoing simultaneously longitudinal and transverse expansion in 113 
involving anisotropic hydrodynamics see 16|). The remarkable agreement between these 
solutions and those of relativistic dissipative fluid dynamics (under the same symmetries) has 
brought great insights about the validity of the hydrodynamic description of the evolution of 
the quark-gluon plasma. However, even in this case the solutions of the AWB equation were 
obtained using iterative numerical methods and it was not known how to obtain analytic 
expressions for the momentum dependence of the single particle distribution function, /, 
and the spatial dependence of its moments. 

nn 

In this paper, we expand on the arguments developed in Ref. p^, 11^ to obtain a new 
fully analytical solution for the single particle distribution function of the AWB equation for 
conformal kinetic systems. The key difference with respect to the exact solutions previously 
derived in [^, 141 involves the global symmetries imposed on the conformal system. The 


symmetry assumptions 


17 


19l | previously employed 


m 


[3,[i4 


14j | were more applicable to the 


matter created in ultracentral relativistic heavy ion collisions while in this work we broaden 
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our focus and consider symmetries more appropriate for conformal systems undergoing three 
dimensional radial expansion, such as the early universe^. We note that the same set of 
symmetries has already been imposed to conformal fluids in 20, l2l| in order to hnd the first 


analytical solutions of second order conformal fluid dynamics. 

The possession of an analytical solution for / allowed us to directly explore here impor¬ 
tant technical aspects in kinetic theory such as the imposition of matching conditions, the 
decomposition of / in its moments in a nontrivial setting as well as its positivity. More im¬ 
portantly, this analytical solution has also revealed a new feature of conformally invariant, 
radially expanding systems described by the AWB equation: the ability to flow as a perfect 
fluid even though the overall dynamics is intrinsically dissipative (e.g., the non-equilibrium 
entropy component is nonzero). In fact, we show that in this solution the energy-momentum 
tensor is exactly that of an ideal fluid at any spacetime point (even though the shear viscosity 
coefficient is nonzero) while the entropy density, computed directly using the full distribu¬ 
tion function, is different than its ideal limit. In this case, this non-equilibrium contribution 
to the entropy density is due to higher order scalar moments (which possess no hydrody- 
namical interpretation) of the Boltzmann equation 22|] that remain out of equilibrium while 
the energy-momentum tensor retains its local equilibrium form. Therefore, in the system 
considered here, slowly moving hydrodynamic degrees of freedom can exhibit true perfect 
fluidity while being totally decoupled from the fast moving, non-hydrodynamical microscopic 
scalar degrees of freedom that lead to entropy production. 

This paper is organized as follows. In the next section we briefly review how AdS 2 <8) S 2 
invariant solutions of fluid dynamics were obtained in Refs. 20|, l21|. In Sec. IIIII we derive 
the main results of this paper and solve the Anderson-Witting-Boltzmann equation for a 
conformal system in AdS 2 ® S 2 geometry. We show in Sec. IIVI how these solutions appear 
from the perspective of the method of moments. We then conclude with a summary of our 
results. Throughout this paper, we use natural units h = c = ks = 1- 


^ An important distinction with respect to the physics of the early universe is that here we still consider an 
underlying flat spacetime. 
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II. RELATIVISTIC HYDRODYNAMICS IN AdS2 O S 2 


We follow 


2l| and consider the ont-of-equilibrium dynamics of a conformal system in 


AdS 2 0 S 2 geometry. This curved geometry is conformally equivalent to 4-dimensional 
Minkowski spacetime (in spherical coordinates), 


ds^ = 


-dt^ + dr^ + r^dn^ 


= — cosh^ p dr"^ -I- dp^ -|- , 


( 1 ) 


where = dO"^ -|- sin^dd0^ is the usual angular piece involving the angles 9 G [0,7r] and 
0 G [0, 27r] while r and p are global AdS 2 coordinates dehned using the Minkowski time, t, 


and 3-dimensional spatial radius, r, in the following way 


2l| 


tanr = 


L2 + r2 - 


cosh p = y/(L2 -K (r -F ty) (L2 -h (r - t ^), 


( 2 ) 


2Lt ’ ^ 2Lr 

with L being the radius of AdS 2 . In this curved space, quantities evolve in r while p plays 
the role of a spatial radial coordinate. In this Weyl rescaled coordinate system the nonzero 
Christoffel symbols are 


= “cosd Sind, = coshp sinhp, r;^ = tanhp. (3) 

The energy-momentum tensor, of a relativistic conformal fluid is usually decomposed 
in terms of the time-like (normalized) local velocity held, as 




Above, we introduced the energy density e = u^UpT^'^, the thermodynamic pressure P (e) = 
e:/3, and the shear stress tensor We further dehned the projection operator 

onto the space orthogonal to A^^" = + u^u'^, and the double, symmetric, traceless 

projection operator A((^ = (AjjA^ -f AjjA^) /2 — A^'^Aa/g/d. Our convention is to dehne 
the huid velocity using the Landau picture, = —eu^, which implies that the energy 

dihusion is always zero. The bulk viscous pressure of a conformal huid is always zero, which 
means that the dissipative processes involving energy and momentum in such systems are 
solely governed by the shear stress tensor. 

The main equations of motion satished by this huid are given by the conservation laws 
of energy-momentum, which we decompose in the following form, 

A^P^T^^ = + A^^a^ InT + A^P^vr^^ = 0, 
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( 4 ) 

( 5 ) 









where -D^ is the general relativistic covariant derivative. The equations above are then 
complemented by the equations of motions for the shear-stress tensor, which, at second 
order in gradients |22l424| . correspond to a relaxation-type equation 

Dxu^ + + higher-order terms , ( 6 ) 

where rj is the shear viscosity and is the shear relaxation time. For a conformal fluid, 
the shear viscosity must be proportional to the entropy, 77 ~ s, while the shear viscosity 
relaxation time must be inversely proportional to the temperature, ~ 1 /T. Above, we 
introduced the shear tensor of the fluid, The brackets () denote the transverse 

and traceless projection of a tensor 

The hydrodynamical solution studied in |2]| was constructed using a static though non- 
uniform local velocity, = (— cosh p, 0, 0, 0) in AdS 2 ( 8 )S 2 space with coordinates (r, p, 9, 0). 
This implies that the system is undergoing a certain type of spherically symmetric radial 
flow in the usual Minkowski coordinates that is equivalent to the conformal soliton flow 


hat was first introduced in 


25l | in the context of the gauge/gravity duality [2^ (see, e.g.. 


2l| for more details about our flow velocity in Minkowski coordinates). With this static 


flow configuration in AdS 2 0 S 2 , the expansion rate of the fluid vanishes, i.e., = 0 , 

and so does the shear tensor, = 0. Thus, Eqs. (Hj) and (jS]) can only be satisfied if the 
temperature and depend solely on the spatial coordinate p, e.g., T (r, p, 9,4)) ^ T (p). 

Moreover, note that in this space is trivial: a quick look at Eq. ([ 6 ]) (and its gener¬ 
alization including terms involving higher order derivatives of the flow) reveals that in this 
problem is identically zero. In fact, since here the flow is static and = 0, = 0, 

and = 7 r^'^(p), in our conformal theory there are no dynamical sources available to induce 
a nontrivial spatial prohle for the shear stress tensor, which must then vanish in all space. If 
nonlinear terms quadratic (or of higher-order) in were present i n , no ntrivial solutions 
of these homogeneous algebraic equations for could be found [ 20 , l 2 l| but those would 
necessarily assume that must be nonzero for any value of p. Therefore, this nontrivial 
branch of solutions is not smoothly connected to the usual hydrodynamic gradient expansion 
for which, in this problem, the first-order Navier-Stokes contribution vanishes. In any case, 
this type of solutions is not going to play a role in our discussion since the nonlinear terms 
in cannot appear in an effective hydrodynamic theory obtained from the Boltzmann 


equation with a linearized collision term 


27l | such as in AWB. Therefore, one can safely set 
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= 0 in the following. Also, since transforms covariantly under Weyl transformations 
24j |. the fact that this quantity vanishes in AdS 2 ^ S 2 implies that it will also vanish in 


Minkowski coordinates. 

In this case, the momentum equation fISl) leads to an equation of motion for the temper¬ 
ature that can be easily solved [21 1 


dpInT = — tanh p T [p) ~ (cosh p) 


-1 


(7) 


The interesting feature of this solution is that it corresponds to the solution of an ideal fluid. 
This happened without making any assumptions about the magnitude of the shear viscosity 
coefficient - it simply appeared as a feature of this highly symmetrical flow conhguration. 
That is, even though the system in principle has a nonzero shear viscosity coefficient, its 
hydrodynamic degrees of freedom cannot dissipate since all gradients are exactly zero in 
AdS 2 0 S 2 (note that dissipation via bulk viscosity is forbidden due to exact conformal 
invariance). In Minkowski space, the temperature evolves in time as it would in a genuine, 
dissipationless fluid. 

In the next sections we investigate the same problem of out-of-equilibrium AdS 2 ( 8 ) S 2 
dynamics from a kinetic theory perspective using the AWB equation. We then clarify which 
non-hydrodynamic degrees of freedom of the microscopic theory are responsible for dissi¬ 
pation in this case and why such degrees of freedom do not couple with the hydrodynamic 
modes. 


III. ANDERSON-WITTING-BOLTZMANN EQUATION 


The on-shell AWB equation in curved spacetime is ISj, [14 1 


P^dpf + T^piPxP^ 


dpi 


P^Up 

Ael 


(/ - /eq) , 


( 8 ) 


m 


where the distribution function / = f{x^,pi) is dehned in a 7-dimensional phase space 
which each point is described by seven coordinates, i.e., the AdS 2 ( 8 )S 2 spacetime coordinates 
= (r, p, 9, 0) and the three spatial covariant momentum components Pi = {Pp,Pe,P(t>)- The 
zeroth component of the momentum is obtained from the on-shell condition for massless par¬ 
ticles PpP^ = 0. Moreover, /eq = exp {p^^Up/T) is the local equilibrium distribution function 
for massless particles with Boltzmann statistics, T is the local temperature, is the local 


6 











velocity of the system, and r^ei is the relaxation time associated with the collision operator. 
Conformal invariance imposes that the relaxation time must be inversely proportional to the 
temperature, Trei = c/T, with c bein g a c onstant that is directly related to the shear viscosity 


27 


29| (thus, the free streaming limit corresponds to 


to entropy density ratio, rj/s = 5c 
c —)■ oo). 

At hrst glance, it may appear that the AWB equation is a linear equation in /. However, 
we note that Eq. ([8]) must be solved simultaneously with the equations of motion for the 
temperature and velocity, Eqs. (j4]) and ([5]). In these, one must also use the dehnition of the 
shear stress tensor of a dilute single component gas, 

cPp 


TT' 


flV _ _ 


(27r)^P^v^ 


/■ 


In the end, one has a coupled set of nonlinear integro-differential equations for /, T, and 
u^. It is commonly very challenging to solve these types of equations even numerically. 


However, as mentioned above, exact solutions o 


obtained using iterative numerical methods 
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this system of equations have been recently 


13 


14j | . For the type of flow and symmetries 


considered in this paper, we demonstrate in the following sections that it is possible to obtain 
analytic solutions of this system of equations. We note that the collisionless limit of a system 
with a flow equivalent to ours in Minkowski space was previously studied in 30| using very 
different techniques than the ones used below. 


A. Analytic Solution 

As mentioned in the previous section when we discussed the fluid dynamical equations, 
the symmetry for the static flow imposes that = (— coshp, 0, 0, 0). Also, for this type of 
static flow / may depend only on the spatial coordinates p, 6*, and 0 (though we shall see 
that / does not depend on this coordinate in the end) and their corresponding momenta. 

Since in the AWB equation the collision term is approximated to be linear in / — /eq, it is 
impossible for terms quadratic or quartic in to appear in the equation of motion for 
at any order in the hydrodynamic series 22, Such terms can only originate from the 


nonlinear terms of the collision operator and, assuming that higher order tensorial moments 
22| initially vanish, the shear stress tensor constructed using the solution / of Eq. (|8]) must 


be zero. Therefore, since the bulk viscous pressure has to be zero due to the underlying 
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conformal invariance, the temperature that enters the AWB equation will satisfy Eq. ([7]) 
with solution 


T{p) = 


Tn 


cosh p ’ 


(9) 


where Tq is a constant. Note that this is not usually the case and in general the temperature 
has to be solved simultaneously with the AWB equation [^, 1^. The fact that the velocity 
prohle is static and the temperature prohle can be solved analytically will be extremely 
useful here since it will allow us to hnd analytical solutions of the AWB equation for this 
system. These solutions for T and serve to considerably simplify the expression for the 
local equilibrium distribution function and the relaxation time, which take the following 
form 


/eq 


Ael 


exp [—cosh p/T (p)] , 
c c 

__ = coshp, 

T [p) To 


( 10 ) 

( 11 ) 


where p'^ = ^Jp^T Pe + {p\l sin^ 6 *) / cosh p. 
The AWB equation then becomes 

Ppdpf - tanhp (^p^ + pi + 

, 1 pI df _ To 1 


pI \ ^ or 

w— + PeOef 
sm 9 J opp 


'pI + pI + 


sin^S 


(/ - /.q) 


( 12 ) 


tan 9 sin^ 9 dpe c cosh p 

where we used Eq. iHD. We note that /eq itself satishes this equation, as is expected for a 
stationary solution (see also the collisionless study of 3^). We also remark that there are no 
terms including df /dp^, which is consistent with spherical symmetry in these coordinates 
and, thus, / does not depend on 0. It is then easy to see that the general solution of this 
equation can be written as a sum of an equilibrium piece and a non-equilibrium part as 
follows: f{p,9-,pp,pe,p^) = /eq feo^{p,9-,pp,pe,p^) where the non-equilibrium piece is 


^{p;Pp,Pn) = J 


^Jpl + pi coshp' 


exp 


.Zk pp_ 

c \Pp\ 


arctan sinh pw 1 -|- 


(13) 


Here, J ( 7 ) is an arbitrary function of its argument 7 and we have dehned the short-hand 
notation = p^T (p^/ sin^ 6 *). By taking c —)■ cxo, one can see that J is actually the solution 
of this equation in the free-streaming limit, J = ‘hfree-streammg- As will be discussed in the 


















following, the functional form of J can be determined by using the matching condition for 
the energy density while requiring that / is positive-dehnite at any point of phase space. As 
far as we are aware, this is the hrst analytical solution of the AWB equation that describes 
a radially expanding system. 


B. Matching condition and positivity 


In kinetic theory it is quite common to dehne the temperature of the system by requiring 
that the energy density of the system is solely determined by its equilibrium value. 


£ = = £eq (T) . 


This condition implies that the following integral must always vanish 


cPp p’" cosh p 


(2i) 


sin6' 


S«i’i{p,0\p„pe,P4) = 0 . 


(14) 


Using the analytic solution derived in the previous section, Eq. (ir3|) . it is possible to reduce 
this integral to a considerably simpler form 


d7 7^jr(7 )exp(-7 ) = 0. (15) 

Now, the condition flT5|) can be used to determine J7’(7). For simplicity, in this work we 
consider a polynomial Ansatz 

5 ^( 7 ) ~ 07 - 1 , (16) 



and one can easily hnd that condition flTSl) is met as long as a = 1/4. Therefore, 


Jii) ~ - 1. (17) 

Note that this function is not positive-dehnite for 7 G [0,4]. However, we still have the 
freedom to hx the overall multiplicative constant. A mandatory physical constraint is that, 
in the end, the distribution function must be a non-negative real-valued function of its 
arguments. In fact, positivity can be obtained as follows. First, note that the sign of the 
exponent in our solution for $, in Eq. (IT^ . is determined by the sign of Pp\ in the limit of 
p —>■ cx), the solution is bounded by exp [— 7rTo/(2c)], when pp > 0, and by exp [7rTo/(2c)], 
when Pp < 0. To make sure that / is positive-dehnite and, at the same time, that hmc_j.o / = 


9 




/eq (i.e., for a vanishing relaxation time one must recover the local equilibrium) we £x the 
overall multiplicative constant to be exp [—To 7 r(l + ^)/{2c)] with .^ > 0 and, thus, 


Jb) = (5 - 1 ) 


exp 


2 c 


■(1 + 0 


(18) 


In principle, other forms of J ( 7 ) may be used in order to achieve the same outcome, which 
would then generate a class of solutions of the AWB equation. In this work, however, we 
limit our discussion to the form fflSj) for 1 /( 7 ). 

It is instructive to study the dependence of / on some of its arguments. For instance, for 

p = 0 

/ 




eq 


p=0 


= 1 + exp 


ttFo 
’ 2 c 


(e + 1 ) 


,2 , .2 , pI 


Jjr\lvi + vt + - 


Arji \l I'D ' I'd ' ' 2 n 

4in V sin 6^ 


(19) 


while for pe = = 0 

f 




eq 


= 1 + exp — 
Pe,P<i>=0 1 2c 


vr(l + ^) + 2 ■p^Gd(pp) 
\Pp\ 


(^bplcoshp- 


where Gd(x) = 2 tan“^ (expx) — 7 r /2 is the Gudermannian function. One can see that these 
expressions are positive-dehnite and that they reduce to the equilibrium distribution in the 
zero mean free path limit c —)■ 0 . 


C. Non-equilibrium entropy 

The local entropy density is computed using the solution for / as follows 

1 f d?p 


s = 


-u^,pP / (In / - 1 ) 


( 21 ) 


(27r)3 y 

It is easy to show that the equilibrium result is Seq = 4 T^(p)/ 7 r^, which is what one would 
expect for an ideal conformal gas with degeneracy factor equal to one. From the form of / 
assumed in this paper, / = /eq (1 + $), one can write the nonequilibrium correction to the 
entropy as 

d?p 


As = s 


Seq — 


The second term can be reduced to 

f d^p 


/eq {(1 + $) In [1 + <F] + $ (In /eq - 1 )} . 






U^,P^ /eq<f> (In /eq - 1) = 


47^2 

T3 

Stt^ 


^ (1 + 7) ^(b) 


( 22 ) 


(23) 


?{(p) exp 


ttTq 
' 2 c 


(1 + 0 
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FIG. 1: (Color online) Relative entropy production As/seq in Eq. (j25h for different values of Tq/c 
(with fixed ^ = 0.01). The solid black line was computed using Tq/c = 10, the dashed red curve is 
for Tq/c = 1, while the dotted-dashed blue curve is for Tq/c = 0.1. 


where 


'H(p) = 2 / dx cosh 


—arctan I 

c V 


f sinh p\ 


X 


(24) 


is a positive-definite function. When going from the hrst line to the second line above, we 
replaced the form of J' ( 7 ) obtained in the previous sections. The full result is 


As 


^eq 


2 -7 
16^ " 


X In <j 1 -F Jipy) exp 
° d7 


dx 


1 + J"( 7 ) exp 

( 


To. _i/sinhp\ 
— tan ' ' 




To _i/sinhp\ 

-tan ' ' 

c 


•1 


16 


7^6'’' 


X In < 1 -F J"( 7 ) exp 


dx 

To 


X J 


1 + J{^) exp 


To 


■ tan 




Wp) 

32 


f sinh p\ 

X J 
ttTo 




exp 


2 c 


(1 + 0 


(25) 


It is easy to see that As(p) is even in p and that As(p) < 0, as expected on physical 
grounds. We show in Fig. [Da plot of As/seq as a function of p for different values of Tq/c. 
For small values of Tq/c one can see that the full entropy density becomes different than the 
equilibrium one (though by a small amount) and that this effect becomes more pronounced 
for large values of p, where it reaches a stationary value that depends on the parameters c 
and Tq. 

Note that As/s^q in Eq. fl25|) does not change under Weyl transformations and, thus, 
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one can find its value in flat spacetime via simple substitution ^(p) = r)). Also, 

since large values of p correspond to large values of t for fixed r (see Eq. ([2])), this quantity 

approaches a constant at large times in flat spacetime. From the result shown in Fig. [H one 

can see that the spatial integral of As in flat spacetime 

128L3T3 r°° r2 As 

^ Jo [L2 + (r + t)2]3/2 _ ^)2]3/2 ~ 

goes to zero when t ^ oo, which indicates that the entropy approaches its equilibrium value 
as time increases. The equation above was obtained using that the equilibrium entropy 
density in flat spacetime is Seq(t,r) = 4:T^(t,r)/{7r'^r^). 



IV. COMPARISON TO THE METHOD OF MOMENTS 


In order to better understand some features of the solution derived in the previous sec¬ 
tions it is convenient to expand $ = (/ — /eq) / fpq in terms of its moments, using irreducible 


tensors and a complete basis of polynomials 22|. The irreducible tensors, 1, 
k''^k''k^\ • • •, are used to expand the angular part of the single-particle distribution func¬ 
tion. They form a complete and orthogonal set, analogously to the spherical harmonics 
jb], and are dehned as k^^^ ...k^^^ = where the transverse, symmetric, and 

traceless projectors are dehned in jb]. Our solution in Eq. flTbll is anisotropic in 

momentum space and hence it possesses both scalar and higher rank moments. For the 
sake of illustration, in this section we focus on the scalar moments of our solution. The 
scalar part of the distribution function is expanded using a set of orthogonal polynomials, 
-^kn ~ 1 where the coefficients aYr were calculated so that 

r dK 


[U^K ) — Or, 


(27) 


22| . Here, we dehned 


(2£-M)!!/ 

using the Gram-Schmidt orthogonalization method as demonstrated in 
dK = d^k/ [(27r)^ k^ and = {—lY/l 2 e^e where, for a nongenerate massless gas of particles, 

(n+1)! T’^+2 
“ (2g + l)!! 27r2 ' 


The irreducible tensors also satisfy orthogonality conditions. 


dK 

7^ 




(2m+ 1)!! 


dK 




V^{2i + 1)\\ 


Fk (u^k^) 


2m 


(28) 
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where Fk is an arbitrary function of u^k^. 

Using this basis, the moment expansion of $ is 

OO OO 

e=0 n=0 

where the moments can be obtained using the orthogonality relations satished by the basis 
elements and are given by 

■■■" = /(2^ (29) 

For the sake of convenience, we defined above 


pW ^ ^ 

' kn £j / j ^mn-^ k 
m=n 

We note that the scalar moments can also be calculated analytically, by replacing Eq. 
(ir3|) into Eq. (I29j). The solution is 

Tl — 2 

©n = r(n + 2) H{p) exp 

where r(n) is the Gamma function. Note that this quantity vanishes for n = 2, as expected, 
from the energy matching condition and that 0o/T^ = ©i/T^ < 0 while > 0 for n > 2. 
For the sake of completeness, in Fig. 2 we plot ©a/T® for different values of Tq/c with 

e = o.oi. 

The actual moment expansion of <F then becomes 


ttTq 
' 2c 


(1 + 0 


(30) 




71=0 


Truncating this expression at n = 2 (note that the matching condition Exes ©2 = 0), we 
obtain something analogous to the 14-moment approximation 23| . 


4 > = pS©o + p£^ 0 i. 


For a gas of nondegenerate massless particles, it is easy to show that 


p(0) ^ ^ 

' kO p2 


3 + ] , pS’ = 


2T ^ 


2^3 


where we used that 


- 1 
“oo ~ -‘-7 


“ii 


1 

*^10 _ _ 9 p 

2T4’ .(0) “ 

“11 
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FIG. 2: (Color online) Normalized scalar moment ©s/T^ for different values of Tq/c (with fixed 
^ = 0.01). The solid black line was computed using Tq/c = 10, the dashed red curve is for Tq/c = 1, 
while the dotted-dashed blue curve is for Tq/c = 0.1. 


In this truncation scheme, the distribution function is then approximated to be 


(l - ^ \/Pp + Te + ( 4 / d) ^ . 

For p = 0, this is exactly the same as our analytical solution, see Eq. ([19]). This shows that 
a hnite number of scalar moments are able to provide a reasonable description of this system 
at least when p = 0. 


$ = -^nip) exp 


ttTq 
’ 2c 


(1 + 0 


V. CONCLUSIONS 

In this paper we derived the hrst analytical solution of the Anderson-Witting-Boltzmann 
equation for a radially expanding system (known as conformal soliton flow) of massless 
particles. We further demonstrated how the matching conditions, commonly used to de- 
hne temperature in kinetic theory, restrict the form of the solution of the single particle 
distribution function. 

The solution we found has some very interesting features. In this system the slowly 
moving hydrodynamic degrees of freedom do not see dissipation, e.g., sound waves propagate 
without any distortion from viscosity. However, faster degrees of freedom are still present 
and they produce a hnite amount of entropy. This may be the hrst example of a kinetic 
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system that does not have a viscous hydrodynamic behavior: between its ideal fluid and 
free-streaming limits, there is no region in space and time where a viscous fluid dynamical 
description is valid. 

This conclusion regarding the perfect fluidity of the conformal soliton flow, studied here 
at weak coupling in the context of kinetic theory, was also found in the case of an infinitely 
coupled AA = 4 Supersymmetric Yang-Mills plasma 25|]. In fact, even though this strongly- 
coupled system has nonzero shear viscosity rj/s = l/( 47 r) Ml, the underlying symmetries of 
the flow together with conformal invariance impose that the energy-momentum tensor of the 
system retains its perfect fluid form. This shows that the exact cancellation of shear viscous 
effects in the energy-momentum tensor discussed here also happens in strongly coupled 
systems. 

For any finite value of c in the relaxation time ffTTl) . our solution for the distribution 
function does not return to local thermal equilibrium even at sufficiently large times. In fact, 
one can see from (j 2 |) that large times (for fixed radius r) correspond to large p’s and, in this 
case, the non-equilibrium contribution given by flTdll and flTSl) does not vanish if c 7 ^ 0 . Thus, 
in our system the effects of the expansion overcome the collision term and the distribution 
function does not relax to its equilibrium form. We note that a similar conclusion was 
found for a different type of rapidly expanding gas in , which went beyond the relaxation 
time approximation and took into account the full nonlinearities of the collision term of the 
relativistic Boltzmann equation. 

The essential approximations made here to And this novel many-body effect were: relax¬ 
ation time approximation, conformal dynamics, and spherical symmetry (implemented via 
the AdS 2 < 8 )S 2 construction). The effects discussed in this paper may appear when describing 
a perfectly radially symmetric and homogeneous droplet of quark-gluon plasma, at very high 
temperatures and vanishing chemical potentials, expanding in vacuum. In this limit, QCD 
is approximately conformal and the flow configuration should resemble the one discussed in 
this paper. 
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